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When treating interactions in quantum dots within a RPA-hke approach, time-reversal symmetry plays an im- 
portant role as higher-order terms - the Cooper series ~ need to be included when this symmetry is present. Here 
we consider model quantum dots in a magnetic field weak enough to leave the dynamics of the dot chaotic, but 
strong enough to break time-reversal symmetry. The ground state spin and addition energy for dots containing 
120 to 200 electrons are found using local spin density functional theory, and we compare the corresponding 
distributions with those derived from an RPA-like treatment of the interactions. The agreement between the two 
approaches is very good, significantly better than for analogous calculations in the presence of time-reversal 
symmetry. This demonstrates that the discrepancies between the two approaches in the time-reversal symmetric 
case indeed originate from the Cooper channel, indicating that these higher-order terms might not be properly 
taken into account in the spin density functional calculations. 
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The impressive progress seen in the recent years in the fab- 
rication and manipulation of semiconductor quantum dotsi 
has been associated with a shift in the perspective with which 
these nanoscopic objects are considered. While the focus used 
to be on the study of their basic properties, most recent ex- 
periments involve them as part of a more complex system in 
which they play a specific role. Some examples include spin 
filtering,^ charge detection,- or the manipulation of quantum 
information.'* 

This new stage of maturity of quantum dot physics makes 
it necessary to go beyond a theoretical understanding of the 
qualitative properties of these systems, and to also develop 
good simulation tools which allow one to predict their behav- 
ior quantitatively. A natural way to proceed is to rely on den- 
sity functional theory.^ '' More precisely, since the spin turns 
out to be an essential degree of freedom, we have in mind a 
spin density functional theory, in the local spin density ap- 
proximation (LSDA), where each spin density is taken as an 
independent variable. The properties of relatively small dots 
have been investigated in great detail^ in this framework, and 
more recently this approach has been used for significantly 
larger dots containing several hundreds of electrons. '^•^ 

In this large dot limit, and assuming chaotic dynamics 
within the dots, it is expected that the fluctuations of ex- 
perimentally relevant quantities, such as the addition en- 
ergy or the ground state spin, should be reliably predicted 
within a Landau Fermi-liquid picture of the Coulomb in- 
teraction and a model of wave functions and eigenener- 
gies based on Random Matrix Theory (RMT) and Random 
Plane Waves (RPW)i'°i"i^^i'^'"* (By "fluctuations" we mean 
the dot-to-dot variation or variation as a function of en- 
ergy typically seen in mesoscopic systems.)' We shall be- 
low refer to these predictions (somewhat inappropriately) as 
the "RPA approach" since the Landau Fermi-liquid picture 
arises from a Random Phase Approximation to the screened 
Coulomb interaction.'" ''' Quite unexpectedly, however, for a 
model chaotic quantum dot the same quantities computed with 
LSDA turned out to differ from these theoretical predictions, 
even at the qualitative leveli^ 



This motivated us to introduce a second-order Strutinsky 
approximation to the LSDA calculations;'^ it was shown to be 
precise enough to serve as a basis for interpreting the LSDA 
results. This analysis indicated that the origin of the discrep- 
ancies between the LSDA and RPA results arises from two 
reasons (on which we shall elaborate below): (i) an effective 
residual interaction which is slightly larger in LSDA; (ii) the 
absence in LSDA of renormalization of the Cooper channel 
associated with higher-order terms in the residual interaction. 
By the Cooper channel we mean the interaction of two elec- 
trons following time-reversed paths for which, then, multiple 
scattering is highly likely.'* 

Because the Cooper channel contribution is non-zero only 
for time-reversal invariant systems, a direct implication of the 
above analysis is that the agreement between LSDA simula- 
tion and RPA prediction should be significantly better in the 
presence of a magnetic field strong enough to break time- 
reversal symmetry. By showing that this prediction actu- 
ally holds, the goal of this paper is to provide an unambigu- 
ous demonstration of the importance of the Cooper channel, 
and as a consequence to indicate a possible limitation of the 
density functional approach to low temperature quantum dot 
properties. 

We shall therefore compare in this paper the distribution of 
addition energies and ground state spins obtained in two dif- 
ferent ways. The first one corresponds to a full fledged LSDA 
calculation for a specific model quantum dot. The second is a 
prediction for a generic chaotic system (i.e. not associated to 
any specific confinement) for which an RPA-Uke treatment of 
the interaction is used. 

The model quantum dot for the LSDA calculations is con- 
fined by the two dimensional potential [r = (a;, y)] 
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and subjected to the uniform magnetic field Bi. (We use 
effective atomic units throughout the paper, which, for GaAs 
quantum dots with effective electron mass m* = 0.067me 
and dielectric constant e = 12.4, corresponds to 10.08 meV 



2 



for energy, and 10.95 nm for length.) Here we use a = 10~^, 
which for particle number in the range [120, 200] gives a 
parameter ~ 1.4 (the ratio between the interaction energy 
and kinetic energy, formally defined as — l/y^Tmoo for 
2D bulk systems). The coupling A between the two oscillators 
is fixed around 0.6 so that the classical motion within the dot 
is in the hard chaos regime. Finally, the parameters b and 
7 eliminate any discrete symmetry. We shall keep b — tt/A 
constant, but, to enhance the statistical significance of the 
various properties studied, we will use the five sets (A, 7) ~ 
{(0.53,0.2), (0.565,0.2), (0.6,0.1), (0.635,0.15), (0.67,0.1)}. 

In the spin density functional description, one considers a 
functional of both spin densities [n"(r), 71'^ (r)] 
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where a — a,j3 correspond to majority and minority spins, 
respectively. The second term is an explicit functional of the 
densities. 



£,,tK,n^]= j dr n{v) U,^t{v) 



fdrdr'^p^+£^,[n'^,nP], (3) 
|r — r'l 

where n{r) = n"{r) + n^{r), Ucxti^) is the exte- 
rior confining potential Eq. Q, and we use the Tanatar- 
Ceperley parametrization of the the exchange-correlation term 

The kinetic energy term TKsin" ,n^], on the other hand, is 
expressed in terms of a set of auxiliary orthonormal functions 
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with A(r) the vector potential generating the magnetic field 

and n'^(r) = Et"il^f(r)l'- 

The LSDA energy for a given value of the total number 
of particles iV = iV" + TV^ and spin S = {N°' - iV^)/2 
is obtained by minimizing the functional Eq. (|3 under the 
constraints / drn°''^{r) = N"'^. The ground state energy for 
a given N is then obtained by picking the spin S with lowest 
energy. The details of the practical implementation of these 
calculations can be found in Ref.ll9l 

Already quite good approximations of the densities n^p{r ) 
are obtained via the Thomas-Fermi approximation (see e.g. 
Fig. 2 in Ref.lish. which amounts to replacing the Kohn-Sham 
kinetic energy Eq. (0} by the explicit density functional 



2N{0) 
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(valid in 2D) with N{0) — m/nh the two-dimensional den- 
sity of states. Once the Thomas-Fermi calculation has been 
performed (see Ref. 20 for the practical implementation), the 
(spin dependent) one-particle effective Hamiltonian 
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FIG. 1: Cumulative density of nearest neighbor level spacing for the 
effective Hamiltonian derived from the Thomas Fermi approximation 
at B = OT (dashed gray) and B = 0.15T (dashed black). The 
solid lines are the corresponding Wigner surmise distributions for 
the orthogonal (gray) and unitary (dark) Gaussian ensembles. 



defined in terms of the Thomas-Fermi self consistent potential 
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gives a good account of the qualitative behavior of the dynam- 
ics within the dot, including self-consistent effects associated 
with the Coulomb interactions. We can therefore use i?TF to 
assess the degree of chaoticity of the dots as well as the ef- 
fectiveness of the time-reversal breaking term. Fig.^displays 
the cumulative density of the nearest-neighbor level spacing 
for the effective Thomas-Fermi Hamiltonian and for two val- 
ues of the magnetic field B = {) and i? = 0.15 T. (The latter 
corresponds to about 7 flux quanta through the area of the dot, 
or a filling factor of about 50 at ^ 1.4.) The perfect agree- 
ment with the random matrix theory prediction for the GOE 
and GUE ensembles, respectively, shows both that the sys- 
tems we consider are indeed chaotic and that B — 0.15 T is 
enough to break time-reversal invariance. We have performed 
additional checks on the statistics of the wavefunctions, and 
in particular on the distribution of the quantities A/y and TVy 
introduced in Eq. (|9|l, which further confirms this point. 

Being confident that the choice of parameters and mag- 
netic field corresponds to the regime we would like to study, 
we can now compute the ground state energy and spin for 
our model quantum dot with the full fledged LSDA self- 
consistent approach. The resulting distributions of addition 
energy and ground state spins, for particle number N in the 
range [120, 200], are shown in Figs.|2]and|3l 

We would like now to compare these distributions with the- 
oretical predictions in which (i) the interactions are treated by 
a RPA-like approach and (ii) the statistical properties of the 
eigenlevels and eigenfunctions are modeled in terms of RMT 
and RPW. The RPA treatment of the interactionsiSiii is essen- 
tially a Landau Fermi-Liquid picture in which the electrons 
are considered to be quasi-particles evolving under an effec- 
tive one-particle Hamiltonian Hcs (which we can think of as 
being Htf) and interacting through a weak short-range resid- 
ual interaction VRpA(r — r')- For broken time-reversal sym- 
metry, VRPA(r — r') can be taken into account by first-order 
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FIG. 2: (Color online) Addition spectra obtained from the LSDA cal- 
culations (dark), compared with theoretical predictions derived from 
RPA/RPW (lighter, blue online). Solid: N even; Dashed: iV odd. 
Top: strength of the RPA interaction obtained by matching the QMC 
value, Js = /o ~ 0.35 (for rs = 1.4). Bottom: better agreement is 
obtained with a slightly larger value of this parameter, Js = 0.46. 



perturbation theory. Up to a smooth component in which we 
are not interested, the energy of the quantum dot for a given 
set of occupation numbers /f = 0, 1 of the orbitals tpi of Hcs 
is therefore the sum of the one-particle energies 
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ip 



and the residual interaction direct plus exchange terms 
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with 



N,. 
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Note that once the interaction potential VFipA(r — r') is de- 
termined, all the relevant quantities are expressed in terms 
of the eigenvalues e, and eigenfunctions ijji of H^tf- The 
RMT/RPW model consists in fixing the eigenvalue fluctua- 
tions according to RMT and those of the eigenfunctions by 



describing them as a superposition of random plane waves, for 
levels within the Thouless energy of the Fermi energy. This 
fixes the fluctuations of Eip and of the Mij and Nij . 

The weak interaction obtained from the RPA treatment'^ — 
is quite naturally the zero frequency and low momentum ap- 
proximation of the screened RPA interaction, which for 2D 
systems takes the form 



^i^PA(r) = 

v^RPA(q) = 



^RPA(q) exp[iq-r] , 
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It is known from quantum Monte Carlo (QMC) results that, in 
the regime of interest, RPA slightly underestimates the inter- 
action strength. Thus we use a boosted interaction 



(11) 



with ^ adjusted so that RPA matches QMC. Note that 
this is completely consistent since the Tanatar-Ceperley 
parametrization of LSDA comes from fitting the QMC results 
as well. The natural quantity to match is the Fermi circle av- 
erage 



Js = (VR,PA)ts — 
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which should be interpreted as the Fermi-Liquid parameter 
/°. For the value — 1.4 considered here, QMC gives^ 
/° = 0.35 while (V[?pA)fe = 0.31. Thus, for our density, the 
natural value is ^ = 1.1; we have used other values as well in 
order to explore trends. 

The distributions of addition energy and ground state spin 
for ^ — 1.1 and 1.5, corresponding to Jg = 0.35 (= /") and 
0.46, are shown in Figs. |2] and |3] We observe that, although 
the lower value of the interaction gives already reasonable 
agreement for the addition energy distribution, a better result 
is obtained, especially for odd, when the interaction is en- 
hanced. In contrast, the spin distribution is accurately repro- 
duced by just taking Js = /q , while the stronger interaction 
predicts slightly too many high spin ground states. 

Despite the apparent contrast here, these results are nicely 
in line with the analysis using the Strutinsky approximation 
in Ref. 15. Indeed, the results there indicate that the main 
source of discrepancy between LSDA and RPA/RPW comes 
from the absence of screening of the Cooper channel in the 
former. Since the Cooper channel does not contribute when 
time-reversal symmetry is broken, the much better agree- 
ment between LSDA and RPA/RPW here, as compared to 
the time-reversal invariant situation considered before,^''^ 
demonstrates the importance of this higher-order interaction 
ejfect in the ground state properties of quantum dots. 

It was also pointed out in Ref. 15 that beyond the effect of 
the Cooper channel, LSDA differed from RPA/RPW in two 
aspects; (i) a slightly larger effective residual interaction, and 
(ii) the existence of spin contamination, i.e. of solutions of the 
Kohn-Sham equations with different wave functions for the 
majority and minority spins. The first of these points explains 
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FIG. 3: (Color online) Spin distribution obtained from the LSDA calculations (left), the RPA/RPW model with = = 0.35 (middle), and 
Js = 0.46 (right). Contrary to the addition spectra, slightly enhancing the RPA interaction does not improve the agreement with the LSDA 
simulation for the spin distribution; rather, it makes it worse. 



why it is necessary to slightly enhance the RPA interaction 
to perfectly reproduce the LSDA addition energy distribution. 
The spin contamination mechanism, however, has a tendency 
to produce slightly smaller ground state spins, and therefore 
offsets the first point when the statistics of spin ground states 
are considered. As a consequence, spin statistics from the 
LSDA calculations nicely agree with the RPA/RPW predic- 
tion without enhancing the RPA interaction. 

To conclude, we have shown that because of the absence 
of the Cooper channel contribution when time-reversal sym- 
metry is broken, the addition energy and ground state spin 
statistics obtained from LSDA calculations for a chaotic quar- 
tic oscillator at weak magnetic field closely follow RPA/RPW 
predictions. This demonstrates in an unambiguous way that 
higher-order interaction effects in the Cooper channel are in- 



deed the main cause of the discrepancies observed between 
LSDA and RPA/RPW in the time-reversal symmetric case. 
The remaining differences can be understood as arising from 
a combination of spin contamination and the slightly larger 
effective interaction that can be derived from LSDA within a 
second-order Strutinsky scheme. Finally, we stress that be- 
yond these differences, both methods predict a great sensitiv- 
ity of the addition energy distribution (at weak magnetic field) 
to the precise value of the exchange parameter Jg. We suggest 
that experimentally measuring these distributions, at a temper- 
ature low enough'^ to resolve the shift in the maxima between 
the odd and even distributions, would provide a sensitive way 
to directly access this quantity. 
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